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We undertake a simultaneous partial wave analysis to proton-proton and neutron-proton scattering data below 
pion production threshold up to LAB energies of 350MeV. We represent the interaction as a sum of delta-shells 
in configuration space below the 3fm and a charge dependent one pion exchange potential above 3fm together 
with electromagnetic effects. We obtain x^ = 2740|;,p + 3835|„y, with a total of A' = 2734|pp +3610|„;, data 
obtained till 2007 and a total number of 46 fitting parameters yielding ;^^/d.o.f = 1.04. Special attention is 
payed to estimate the errors of the phenomenological interaction as well as the derived effects on the phase- 
shifts and scattering amplitudes. 
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The NN interaction plays a central role in Nuclear Physics 
(see e.g. Ill, 12] and references therein). The standard proce- 
dure to constrain the interaction uses a partial wave analy- 
sis (PWA) of the proton-proton (pp) and neutron-proton (np) 
scattering data below pion production threshold {3}. The 
best PWA uses a large body of NN scattering data giving 
X^/d.o.f < 1 after discarding about 20% of 3(7 inconsistent 
data 1 4]. This fit incorporates charge dependence (CD) for the 
One Pion Exchange (OPE) potential as well as electromag- 
netic, vacuum polarization and relativistic effects, the latter 
being key ingredients to this accurate success. The analysis 
was more conveniently carried out using an energy dependent 
potential for the shoit range part. Later on energy indepen- 
dent high quality potentials were designed with almost identi- 
cal x^/d-o.f ^ 1 for the gradually increasing database DSHil. 
While any of these potentials provides individually satisfac- 
tory fits to the available experimental data, an error analysis 
would add a means of estimating quantitatively the impact of 
NN-scattering uncertainties in Nuclear Structure calculations. 
In the present work we provide a new high-quality potential 
implementing an analysis of its parameter uncertainties using 
the standard method of inverting the covariance matrix [9]. 

The complete on-shell NN scattering amplitude contains 
five independent complex quantities, which we choose for 
definiteness as the Wolfenstein parameters H 

M(k/,k/)=fl + m((7i,n)((72,n) + (g-/!)(c7i,m)(CT2,m) 

+ (g + /z)(ai,l)((72,l)+c(c7i + (72,n), (1) 

where a,m,g,h,c depend on energy and angle. To determine 
these parameters and their uncertainties we find that a conve- 
nient representation to sample the short distance contributions 
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to the NN interaction can be written as a sum of delta-shells 
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«=l 



^V,-.„r,-5(r-r,) 
(=1 



VoPEW+VemW e{r-rc) 



(2) 



where 0„ are the set of operators in the AVI 8 basis fl, 
r, < re are a discrete set of A^-radii and V,.„ are unknown co- 
efficients to be determined from data. The r > re piece con- 
tains a CD OPE and electromagnetic (EM) corrections which 
is kept fixed throughout. The solution of the corresponding 
Schrodinger equation in the (coupled) partial waves ^■^+'Ly 
for r < re is straightforward since the potential reads 



VfHr) = 



1 



2Ala/3 t 



^(A,)/,^5(r-r,) 



r < rr 



(3) 



with pLa^ = MaMp/{Ma +Mp ) the reduced mass with a,fi = 
n,p. Here, (A;)^"^, are related to the V,.„ coefficients by linear 
transformation at each discrete point r,. Thus, for any r, < r < 
r,+i we have free paiticle solutions and log-derivatives are dis- 
continuous at the r,-points so that one generates an accumu- 
lated S-matrix at any sampling point providing a discrete and 
purely algebraic version of Calogero's variable phase equa- 
tion uE\. 

This form of potential effectively implements a coarse 
graining of the interaction, first proposed 40 years ago by 
Aviles 1 17]. We have found that the representation (|3]l is ex- 
tremely convenient and computationally cheap for our PWA. 
The low energy expansion of the discrete variable phase equa- 
tions was used already in Ref. lUSll to determine threshold pa- 
rameters in all partial waves. The relation to the well-known 
Nyquist theorem of samplin g a signal with a given bandwidth 
has been discussed in Ref. lll9ll . Some of the advantages of 
directly using this simple potential for Nuclear Structure cal- 
culations have also been analyzed ll20ll . 

The fact that we are coarse graining the interaction enables 
to encode efficiently all effects operating below the finest 



TABLE I. Fitting delta-shell parameters (^n )/'J/ (in fm ) with their errors for all states in the JS channel. We take N = 5 equidistant points with 
Ar = 0.6fm. — indicates that the corresponding {X„)^, = 0. In the first line we provide the central component of the delta shells corresponding 
to the EM effects below re = 3fm. These parameters remain fixed within the fitting process. 



Wave 


A, 
(ri = 0.6fm) 


A2 
(r2= 1.2fm) 


(r3 = 1.8fm) 


ir4 = 2.4fm) 


(r^ = 3.0fm) 


Vc[pp]em 


0.02091441 


0.01816750 


0.00952244 


0.01052224 


0.00263887 



'5o[pp] 


1.32(2) 


-0.720(4) 


-0.189(2) 


- 


-0.0215(3) 


'5o[np] 


1.37(7) 


-0.77(2) 


-0.17(1) 


- 


-0.025(1) 


'Po 


- 


1.11(2) 


-0.371(7) 


-0.041(3) 


-0.030(1) 


IPI 


- 


1.27(2) 


- 


0.069(2) 


- 


'h 


- 


1.379(5) 


- 


0.0569(5) 


- 


'Sr 


1.52(6) 


-0.44(1) 


- 


-0.061(1) 


- 


£1 


- 


-1.64(1) 


-0.38(2) 


-0.206(8) 


-0.024(3) 


'D, 


- 


- 


0.38(1) 


0.081(9) 


0.016(3) 


'D2 


- 


-0.23(1) 


-0.200(3) 


- 


-0.0198(2) 


'D2 


- 


-0.92(4) 


-0.23(2) 


-0.210(8) 


-0.019(2) 


'Pi 


- 


-0.485(1) 


- 


-0.0241(6) 


-0.0061(3) 


£2 


- 


0.25(2) 


0.212(4) 


0.040(2) 


0.0162(5) 


'F2 


- 


3.47(6) 


-0.234(4) 


- 


-0.0139(6) 


1F3 


- 


- 


0.17(2) 


0.071(7) 


- 


'D, 


- 


0.55(2) 


- 


- 


- 



TABLE n. Deuteron static properties compared with empirical values and high-quality potentials calculations 



Delta Shell Empirical[ 10-15] Nijml[5] Nijm II [5J 



Reid93 [5] 



AVI 8 [6] 



CD-Bonn [7] 



£rf(MeV) 


Input 


2.224575(9) 


Input 


Input 


Input 


Input 


Input 


V 


0.02493(8) 


0.0256(5) 


0.02534 


0.02521 


0.02514 


0.0250 


0.0256 


A5(fmi/2) 


0.8829(4) 


0.8781(44) 


0.8841 


0.8845 


0.8853 


0.8850 


0.8846 


r„,(fm) 


1.9645(9) 


1.953(3) 


1.9666 


1.9675 


1.9686 


1.967 


1.966 


Qoifrr?) 


0.2679(9) 


0.2859(3) 


0.2719 


0.2707 


0.2703 


0.270 


0.270 


Pd 


5.62(5) 


5.67(4) 


5.664 


5.635 


5.699 


5.76 


4.85 


(r-')(fm-i) 


0.4540(5) 






0.4502 


0.4515 







resolution Ar which corresponds to the shortest de Broglie 
wavelength corresponding to the pion production threshold, 
Amin ^ ^/VinjiMN ^ 0.55fm, so that a maximal number of 
points A^ = rc/Ar ^ 5 (for r^ = 3fm) should be needed. In 
practice, we expect the number of sampling points to decrease 
with angular momentum as the centrifugal barrier make ir- 
relevant those points r; < (/ + 1/2)//? below the relevant im- 
pact parameter, so that the total number of integration points 
and hence fitting strengths V,,, will be limited and smaller than 
N = 5. 

The previous discretization of the potential is just a way 
to numerically solve Schrodinger equation for any given po- 
tential where one replaces V{r) — > V{r) — L,- V(r,)r/5(r — r,), 
but the number of delta-shells may be quite large for fixed 
strengths y(r,). For instance, for the ^Sq wave and for the 
AVI 8 1 6] potential one needs N=600 delta-shells to reproduce 
the phase-shift with sufficient accuracy (below 10^^ degrees) 
but just A^ = 5 if one uses V,- as fitting parameters to the same 



phase shift 12C 

The EM part of the NN potential gives a contribution to the 
scattering amplitude that must be taken into account properly 
in order to correctly calculate the different observables. Each 
term of the electromagnetic potential in the pp and np chan- 
nels needs to be treated differently to obtain the corresponding 
parts of the total EM amplitude. The expresions for the con- 
tributions coming from the pp one photon exchange poten- 
tial Vci , and the corresponding relativistic correction Vci, are 
well known and can be found in |4]. To calculate the con- 
tribution of the the vacuum polarization term Vvp we used 
the a ppro ximation to the amplitude given in ||2lll . Finally, 
Ref. 12211 details the treatment of the magnetic moment in- 
teraction Vmm for both pp and np channels and the necesary 
corrections to the nuclear amplitude coming from the electro- 
magnetic phaseshifts. 

Of course, once we admit that the interaction below re 
is unknown there is no gain in directly extending the well- 
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FIG. 1. Deuteron Forni Factors with error bands by propagating the uncertainties. Note the tiny width of the bands. 



known charge-dependent OPE tail for r < re- Unlike the 
purely strong piece of the NN potential the electromagnetic 
contributions are known with much higher accuracy and to 
shorter distances (see e.g. Ref. H) so that one might extend 
Vem{f) below re adding a continuous contribution on top of 
the delta-shells, so that the advantage of having a few points 
in the region r < re would be lost. To improve on this we 
coarse-grain the EM interaction up to the pion production 
threshold. Thus, we look for a discrete representation on 
the grid of the purely EM contribution Vem('"), i-e. we take 
Vemir) = LnVfndir - r.) + e{r - r,)V,„,{r), where the Vf 
are determined by reproducing the purely EM scattering am- 
plitude to high-precision and are not changed in the fitting 
process. The result using the EM potential of Ref. la] just 
turns out to involve the Coulomb contribution in the central 
channel and is given in the first line of Table J] As expected 
from Nyquist sampling theorem, we need at most N — 5 sam- 
pling points which for simplicity are taken to be equidistant 
with Ar = 0.6fm between the origin and r^ = 3fm to coarse 
grain the EM interaction below r < re- Thus we should have 
vf = V"'' +Vf if charge symmetry was exact in strong inter- 
actions for r < re, although some corrections are expected as 
documented below. 

In our fitting procedure we coarse grain the unknown short 
range part of the interaction from the scattering data. We use 
a large database compiling proton-proton and neutron-proton 
scattering data obtained till 2007 ll23l - l25ll and carry out at any 
rate a simultaneous pp and np fit. Unfortunately, some groups 
of these data have a common but unknown normalization. We 
thus use the standard floating by including an additional con- 
tribution to the ;^^ as explained in detail in Ref. |i8J. We also 
apply the Nijmegen PWA |4] 3(J-criterion to reject possible 
outliers from the main fit with a 3(7-confidence level, a strat- 
egy reducing the minimal x^ but also enlarging the uncertain- 
ties. Initially we consider A^ = 2900|pp + 5034|„p data and 

get xiin - 3375|pp + 9217|„p yielding zVd-O-f- = 1-59. Ap- 
plying the 3ff-rejection we discard 1590 data. After re-fitting 
the remaining A^ — 2734-\pp + 3610|„p data we finally obtain 
X^i„ = 2140\pp + 3835|„p with a total number of 46 fitting 
parameters so that ;|^^/d.o.f. = 1.04. 

While the linear relations of the (A,)^'^, and V, „ parameters 
are straightforward, limiting the number of operators 0„ re- 
duces the number of independent partial waves. The fitting 



parameters (A„)^'J, entering the delta-shell potentials as inde- 
pendent variables, Eq. (|3]i, are listed in Table I] with their de- 
duced uncertainties. All other partial waves are consistently 
obtained from those. Our final results allow to fix the same pp 
and np potential parameters with the exception of the central 
components of the potential as it is usually the case in all joint 
ppH-np analyses carried out so far |0-|7|. 

We find that introducing more points or equivalently reduc- 
ing Ar generates unnecessary coiTelations and does not im- 
prove the fit. Also, lowering the value of r^- below 3fm, re- 
quires overlapping the short-distance potential, Eq. Q, with 
the OPE plus em coiTections. We find that independent fits 
to pp and np, while reducing each of the j^-values, drive 
the minimum to incompatible parameters and erroneous np 
phases in isovector channels. Actually, the pp data constrain 
these channels most efficiently and in a first step pp-fits where 
caiTied out to find suitable starting parameters for the corre- 
sponding np-phases. Quite generally, we have checked that 
the minimum is robust by proposing several starting solutions. 

As a numerical check of our construction of the amplitudes 
we reproduced the Wolfenstein parameters for the Reid93 and 
Nijmll potentials to high accuracy using A^ ~ 12000 delta- 
shells grid points, which ensures coiTectness of the strong 
contributions. As a further check of our implementation of 
the long-range EM effects along the lines of Refs. |4, 2ll 12211 
we have also computed the ;if^/d.o.f. for Reid93, Nijmll and 
AV18 potentials (fitted to data prior to 1993) which glob- 
ally and bin-wise are reasonably well reproduced when our 
database (coinciding with the one of Ref. ISJ for np) includes 
only data prior to 1993. 

We determine the deuteron properties by solving the bound 
state problem in the ^Si —^ D\ channel using the correspond- 
ing parameters listed in Tablel] The predictions are presented 
in table HI] where our quoted etTors are obtained from prop- 
agating those of Table |T] by using the full covariance matrix 
among fitting parameters. The comparison with experimental 
values or high quality potentials where the deuteron binding 
energy is used as an input is satisfactory [4-8]. 

The outcoming and tiny error bands for the Deuteron form 
factors (see e.g. |26]) are depicted in Fig. [T] The rather small 
discrepancy between calculated and experimental values is 
statistically significant and might be resolved by the inclusion 
of Meson Exchange Currents. In Fig.|2]we show the active pp 
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FIG. 2. np and pp phase shifts and their errors (solid band) corresponding to independent operator combinations corresponding to the fitted 
potential. The points with eii'or bars represent the mean value and standard deviation of the compilation of the PWA L4J and the six high quality 
potentials l^^. 
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FIG. 3. np (left) and pp (right) Wolfenstein parameters (in fm) as a function of the CM angle (in degrees) and for iiLAB = 200MeV. We compare 
our fit (black band) with a compilation (blue band) of the PWA and six high quality potentials I4i,i6r-^ which provided a pj /d.o.f < 1. 



and np phases in the fit with their propagated errors and com- 
pare them with a compilation of the PWA and six high quality 
potentials [4,^-^ which provided a ;i^^/d.o.f < 1. Likewise, 
in Fig.|3]we also show a similar comparison for the pp and np 
Wolfenstein parameters. 

Finally, as the previous analyses [l4j-t8i|] and the present paper 
show the form of the potential is not unique providing a source 
of systematic errors. A first step along these lines has been 
undertaken in Ref. ll27ll . Thus, the uncertainties will generally 
be larger than those of purely statistical nature estimated here. 

To summarize, we have determined a high-quality proton- 



proton and neutron-proton interaction from a simultaneous 
fit to scattering data and the deuteron binding energy with 
X^/A.o.f. = 1.04. Our short range potential consists of a 
few delta-shells for the lowest partial waves. In addition, 
charge-dependent electromagnetic interactions and one pion 
exchange are implemented. We provide error estimates on our 
fitting parameters. Futher details will be presented elsewhere. 

We warmly thank Franz Gross for useful communications 
and providing data files. We also thank R. Schiavilla and R. 
Machleidt for communications. 
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